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Inverse Geometric Approach to the Simulation of the Circular Growth. 
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We demonstrate the power of the genetic algorithms to construct the cellular automata model 
simulating the growth of 2-dimensional close-to-circular clusters revealing the desired properties, 
such as the growth rate and, at the same time, the fractal behavior of their contours. The possible 
application of the approach in the field of tumor modeling is outlined. 
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Introduction 

The surface phenomena, such as the growth, rough- 
ening and penetration into the surrounding medium, as 
well as their modeling and efficient control, represent 
long standing research topic in many branches of sci- 
ence and technology. The special case is the close-to- 
circular growth observed very often in biological systems, 
like bacterial colonies and tumors. The growth mod- 
els are usually classified into two groups, (1) continuum, 
formulated through differential equations (DE), and (2) 
discrete lattice models, most often represented by cel- 
lular automata (CA) agent-based 0, Q, and Monte 
Carlo inspired models. The CA models are typically con- 
structed by a priori (or microscopically reasoned) setting 
up the transition rules ( CA rules) generating the desired 
behavior. However, this ad hoc construction of the CA 
rules is, due to high nonlincarity and multidimcnsional- 
ity of the CA models, often not appropriate in the case 
of more complicated behaviors. Below we propose the 
application of genetic algorithms (GA) Q to retrieve the 
CA rules generating the growth of the clusters revealing 
the desired features, such as the growth rate and fractal 



dimension of the contour. In the here presented inverse 
geometric approach we quantify the differences between 
the properties of a CA cluster growth and the fixed de- 
sired properties (i. e. the quality of the corresponding 
CA rules) by the definition of the objective function, and 
use the optimization technique to set the CA rules which 
minimize the objective function value. We demonstrate 
that the approach is, in general, efficient enough to re- 
trieve the CA rules generating the growth of the CA clus- 
ter with the desired properties. At the end we outline its 
possible contribution to the cancer research. 



Methods 

Cellular automata |5[ were originally introduced by 
John von Neumann as a possible idealization of biological 
systems. In the simple case they consist of a 2D lattice 
of sites i,j = 0,..., Ld, in the states s\f <E {0, 1}, where 
t = 0, . . . , r, is the discrete time and Ld the size of the 
lattice. During the r time steps they evolve obeying the 
set of local transition rules (CA rules) a, formally written 
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which defines the CA rules a as the mapping 

a: {0,f}x{0,I}x...x{0,f}^{0,f}. (2) 



The two basic update strategies are typically adopted: 
i) synchronous, when the t + 1 states ([TJ of all the CA 
sites are computed one by one in advance and the up- 
dates themselves are performed afterwards; and ii) asyn- 
chronous, when a randomly chosen site is immediately 
updated. The CA evolution is represented by the re- 
spective point in 2 9 -dimensional binary space enabling, 
in principle, the enormous number of 2 512 possible global 
behaviors, predestining the CA to be an efficient simula- 
tion and modeling tool [||. Inherent nonlinearity of the 
CA models is, however, a double-edged sword. On the 
one hand, it enables to model a broad variety of behav- 
iors, from trivial to complex, on the other hand it results 
in difficulty with retrieving the transition rules generat- 
ing the desired global behavior. 

No well-established universal technique exists to solve 
the problem, and, despite sporadic promising applica- 
tions of the genetic algorithms to solve the task @,H, Q, 
one typically implements the CA by intuitive (ad hoc or 
microscopically reasoned) construction of the transition 
rules, or reduces the dimensionality of the CA search 
space by incorporating a symmetry into the transition 
rules. Well known are totalistic rules, when only the num- 
ber of the adjacent sites in the state 1 (or 0) affects the 
future state of the respective site. Despite the crucial 
simplification, they provide very complex behavior (e. g. 
Conway's Game of Life). 

Genetic Algorithms [J] are general-purpose search and 
optimization techniques based on the analogy with Dar- 
winian evolution of biological species. In this respect, 
the evolution of a population of individuals is viewed 
as a search for an optimum (in general sense) genotype. 
The feedback comes as the interaction of the resulting 
phenotype with environment. Formalizing the basic evo- 
lutionary mechanisms, such as mutations, crossing-over 
and survival of the fittest, the fundamental theorem of 
GA was derived (schema theorem) which shows that the 
evolution is actually driven by multiplying and combin- 
ing good (quantified by an appropriate objective func- 
tion), correlations of traits (also called schemata, or hy- 
perplanes). The remarkable feature resulting from the 
schema theorem is the implicit parallelism stating that 
by evaluating a (large enough) population of individuals, 
one, at the same time, obtains the information about the 
quality of many more correlations of the traits. 

In the following we use the GA optimization to retrieve 
the CA rules producing the growth of the 2D CA cluster 
by the required rate as well as fractal behavior of the 
contour, which is here understood as the set of all the 
CA sites in the state I with at least one neighbor in the 
state 0. 

At the beginning, the chain of concentric circles (learn- 
ing patterns) with randomly deformed close-to-circular 

contours pf) 6 {0, 1}, i, j = 0, . . . , Ljj, for t = 0, . . . , r, 



were generated accordingly to 
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or 1 drawn with probability 
for \^ Ji 2 + j 2 - RW\ < f, 
for y/i 2 + j 2 < i?W - f , 



with linearly increasing tumor radius 
rW = R<f» + At, 



(3) 



(4) 



where R^ is the initial radius of the CA cluster and A 
(< 1) the desired growth rate constant. 

The first objective of the optimization task solved by 
the GA was to retrieve the CA rules © providing the 
growth from the initial pattern {s^ }, generated by the 

same relation as {p^ } ^ , through the sequence of {s fj } 
([T]). t = 1. . . . .t. with the minimum mismatch from the 
learning patterns {Pif} in the respective t, quantified by 
the objective function 
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(5) 



where S is the standard Kronecker delta symbol, wq is the 
weight factor, in our case wq = 2. The above expression 
of the objective function |5]) reflects the programming 
issues. The larger overlap of {s^} with {p^ } enhances 
the denominator of Eq. ([5]) , the prefactor p^ in the term 
PijS PijSij reduces the computational overhead by ignoring 
the calcul S PtjStj for p y = 0. 

The second objective to the desired growth relates to 
the geometric properties of the contour. Broadly ac- 
cepted invariant measure expressing the contour irreg- 
ularity is the fractal dimension, Dp. Using the box- 
counting method it can be calculated as 



Df = lim 



log M B (e) 



o log(l/e) 



(6) 



where Mb(e) is the minimum number of boxes of size e 
required to cover the contour. Here, it has been deter- 
mined as the slope in the log- log plot of -Ms(e) against 
1/e using the standard linear fit. To obtain the CA rules 
generating the cluster with the desired fractal dimension, 
D-p, the objective function ([5]) has been multiplied by the 
factor 



f 2 (a) = l + Wl (D T F -D ¥ ) 2 , 



(7) 



where Dp is the fractal dimension of the cluster kept after 
the r steps with the CA rules a, w\ being the weight. 
Finally, the objective function is written 

m=Ma)Ma), (8) 

and the optimum CA rules, a*, meet the requirement 

/(0 = min/(a). (9) 
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Results 

If not stated otherwise, all the below CA runs started 
from the pattern {jpy } ©, with the radius = 5. 
The size of the 2D CA lattice Lr> varied form 120 to 
200 (reflecting the growth constant), with the periodic 
boundaries, and the length of the CA evolution r = 200 
steps. The asynchronous update has been adopted. The 
canonical GA search has been applied to retrieve the set 
of the optimum CA rules, a* , which gives minimum ob- 
jective function values (JHJ or [5J respectively). The size of 
the population was kept constant (1000 individui), the 
probability of bit-flip mutation 0.001, and the crossing 
over probability 1. The ranking selection scheme was 
applied. The length of the optimization was 3000 gen- 
erations. Figs. Q] to 2] show the capability of the above 




FIG. 1: The comparison of the typical contours created by 
the CA rules minimizing the criterion (J3J) for the three dif- 
ferent required rates (mediated by the rate constant A = 
0.2,0.3,0.4 in the Eq. Q) (up, bold line) against the pre- 
scribed horizons after t — 0, 50, 100, 150, 200 CA steps. Below 
are the respective average contours calculated from 100 repe- 
titions of the growths. The criterion Q has not been applied 
(wi =0). 

approach to fulfill the required objectives. In Fig. [TJwe 
present the comparison of the typical growths created by 
the CA rules minimizing the criterion ([5]) with the respec- 
tive prescribed rates (|4]) as well as the respective average 
contours calculated from 100 repetitions for each of the 
growths. One can see that in this case the CA rules pre- 
serve circular symmetry for all the tested growth rates. 
Fig. [2] shows the growing contours created by the CA 
rules minimizing both the criteria, the desired rate and 
fractal dimension of the contour, as implemented in 
The respective average contours (Fig. [3]) shows slight de- 
viations from the circular symmetry. The above results 
demonstrate that the GA optimization procedure can be 
used to retrieve the CA rules generating the close-to- 
circular growths by desired rates and, at the same time, 
fractal dimensions of their the contours (Fig. [4]). 




FIG. 2: The comparison of the typical contours created by 
the CA rules minimizing the criterion ([8} for the four different 
required fractal dimensions (Df — 1.1, 1.2, 1.3, 1.4, denoted as 
A, B, C, D) (up, bold line) with the prescribed horizons after 
t = 0, 50, 100, 150, 200 CA steps. In Fig. |3] are the respective 
average contours calculated from 100 repetitions. The crite- 
rion (0 has been applied with the weight w\ = 10. The rate 
constant A in Eq. (|4]) has been kept 0.3 in all the cases. 




FIG. 3: The average contours calculated from 100 repetitions 
for each of the A-D growth in Fig. [2] The deviation from the 
circular symmetry is obvious, especially if the lower fractal 
dimensions are required during the optimization. 
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FIG. 4: Comparison of the average fractal dimensions of the 
contours, Dp, produced by the CA rules optimized to mini- 
mize the objective function (|8]) against the respective desired 
fractal dimension (Df = 1.1, 1.2, 1.3, 1.4, and 1.5). The cir- 
cles show the fractal dimensions for the respective CA rules 
and time step. In all the cases A = 0.3 and wi — 10 (Eqs. [4] 
and [7]). 

Application to tumor growth modeling 

Here we present the application of the above approach 
in the field of tumor modeling, enabling to study the 
tumor growth on the spatio-temporal scales where the 
microscopic details have vanished and eventual univer- 
sal features of the growth emerge. To understand the 
universal properties of the growth, physicists have devel- 
oped the fractal and scaling analysis [lfj , enabling them 
to classify many growth processes of seemingly different 
nature into the universality classes fill ] . Led by similar 
motivations, the fractal and scaling methodology seems 
to be potentially promising approach to study various 
aspects of carcinogenesis. The morphology of growing 
tumor surface is the product of the interactions of many 
biochemical agents, physical processes and geometrical 
constraints. Pathologists have used the morphology and 
shapes of neoplasms to distinguish between the malig- 
nant and benign forms as a rule of thumb. Recently, a 
few papers [lj[ EH EH EH EH have indicated that the 
fractal dimension of the tumor contour, morphology or 
vascular networks can be applied to distinguish between 
the normal and malignant tissue. In vitro grown tumors 
usually form 2 (or 3) dimensional close-to-spherical ag- 
gregations, called multicellular tumor spheroids (MTS) 
|17| . They can be viewed as the prototype of some real 
types of human tumors, such as melanomas or gliomas. 
Allowing the strict control of nutritional and mechanical 
conditions, the MTS are excellent experimental objects 
to test the validity of the models of tumor growth. 

The below outlined application of the inverse geomet- 
ric approach to tumor growth modeling is based on the 
morphometric MTS data by Bru et al. [H, Ell ; namely 



therein suggested linear increase of the MTS mean ra- 
dius, proliferation activity located at the tumor border 
and the range of the fractal dimensions of the in vitro 
grown cell colony contours between 1.1 and 1.4, such 
as HT-29 cell line (tumor type colon adenocarcinoma, 
Dp = 1-13), HeLa (cervix carcinoma, 1.25), Saos-2 (os- 
teosarcoma, 1.34), etc. Using the optimized CA models, 
we reproduce and interpret some of the features, like the 
growth rate, response of the tumor to chemotherapy and 
scaling properties of the contour. We emphasize that in 
the below application the CA site does not correspond to 
the only biological cell. We interpret it as a 2D constant 
size window occupied by the above (or under) critical 
number of tumor cells, corresponding to 2 states - 1 and 
0, respectively. Switching the state of the CA site from 
to 1 (as well as from 1 to 0) during the simulation sim- 
ply means that the number of tumor cells in that site 
has just crossed the critical value. The mechanism itself 
stays out of the scope of this phenomenological approach, 
nevertheless it encompasses the mean effect of proliferat- 
ing activity, cells migration, apoptosis, immune response, 
etc. 



Radius increase 

The simplest mathematical models of the tumor 
gro wth assume an exponential increase of the tumor mass 
[20j | (which is however plausible for the early stage of 
carcinogenesis) . Bru et al. [H, Eil have recently experi- 
mentally shown that the MTS mass, actually grows 
linearly increasing its mean radius, 

= 7r(i?(°) + At) 2 , (10) 

where and A retains their meanings 

Fig. [5] shows the log-log time dependence of the radius 
of the two growing CA clusters applying the different 
rules, optimized to produce the clusters with the frac- 
tal dimension of their contours, Dp, equal 1.1, and 1.5, 
respectively. Their slopes differ negligible, nevertheless 
both are slightly steeper than the slope of the idealized 
growth accordingly to Eq. (jTUJ) . 

Intuitively, the assumption of the linear increase of the 
radius is consistent with the observation that prolifer- 
ation activity is located at the tumor contour. Never- 
theless, the results (Fig. [5]) show slight deviations from 
the linear increase of the radius, which we attribute to 
longer than ideally circular contour of the cluster, giving 
higher chance to proliferate. On the other hand, the frac- 
tal dimension of the contour, Dp, corresponding to the 
roughness of the cluster contour, does not influence the 
growth rate in a simple way. We have observed that the 
CA cluster enforced to grow with the fractal dimension 
of the contour close to 1 very often "freeze" repeating 
some simple pattern which leads to significant deviations 
from circularity f2lj . As the deviations from the circular- 
ity of the tumor contour is the obvious experimental fact 
as well, we find the hypothesis of linear increase of the 
cluster mean radius valid for the spatio-temporal scales 



1 1.5 2 2.5 3 

log (t) 

FIG. 5: Comparison of the growing CA cluster radius created 
by the two CA rules producing the contours with the different 
fractal dimensions Dp = 1.1 and 1.5, respectively, against the 
required case (Eq. I10[) . In both the cases the rate constant 
,4 = 0.3. 

where the deviations of the contour points from circular- 
ity is negligible in comparison to the mean radius. 

Simulation of chemotherapy 

To outline the possible therapeutic perspective of the 
above inverse geometric approach we present the simula- 
tion of the impact of chemotherapy to MTS growth. Most 
chemotherapeutic drugs effectively damage fast prolifer- 
ating cells (such as tumor cells). Here, the chemother- 
apy is included into the model in a similar way as in DE 
models |22j as the term inhibiting the production of tu- 
mor cells, regardless of the microscopic background of the 
process. For our purposes, the chemotherapy is viewed as 
any mechanism (additional to those already imprinted on 
the CA rules by the optimization procedure) converting 
the CA site from 1 (i. e. occupied by the above criti- 
cal number of tumor cells) to (occupied by the under 
critical number of tumor cells) with non-zero probabil- 
ity, represented by the free parameter, called efficiency 
of chemotherapy and denoted as c. In our CA simulations 
the above mechanism was realized simply by switching a 
randomly chosen CA site (of the contour only - by defi- 
nition in the state 1) to with the probability c, instead 
of applying the respective CA rule. 

Fig. [5] shows the typical impact of the simulated 
chemotherapy on the growth of the CA clusters. One 
CA cluster was grown with the transition rules which had 
been optimized to generate smoother contours (quanti- 
fied by its fractal dimension, Dp = 1.1 , upper line), 
one CA cluster grown with the rules optimized to gen- 
erate more rough contours (Dp = 1.4, bottom line). In 
the left column we compare the sizes and shapes of both 



FIG. 6: Simulation of the impact of chemotherapy to the CA 
cluster after 200 time steps of the growth. The left column 
shows the CA clusters grown accordingly to the CA rules op- 
timized to generate the clusters with the fractal dimension 
of their contours, Dp, equal to 1.1 (up), and 1.4 (bottom), 
respectively, without chemotherapy. The right column de- 
picts the cases received with the respective rules, with the 
chemotherapy added (c = 0.15). 

the clusters after 200 time steps grown in the absence 
of chemotherapy (c = 0), against the respective cases 
with the chemotherapy added (c = 0.15, right column). 
Fig. [7] shows the dependence of the CA cluster growth 
on the efficiency of the chemotherapy for the two sets of 
transition rules (generating the contours with Dp = 1.1, 
and 1.4, respectively) in general. The simulation results 
clearly demonstrate that the impact of chemotherapy is 
more inhibiting to the CA clusters grown with smoother 
contours than those growing with more rough contours. 




c (efficiency of chemotherapy) 

FIG. 7: Dependence of the impact of chemotherapy to the 
evolution of the CA clusters generated by the two groups of 
the CA rules. The first set consists of 6 different CA rules 
producing (in the absence of chemotherapy) the clusters with 
the fractal dimension of their contours 1.1, the latter of 6 CA 
rules producing the clusters with the fractal dimension 1.4. 
The CA rules from the latter group produce the CA cluster 
growth less sensitive to the chemotherapy. 
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Hopefully, the above demonstrated results could moti- 
vate researchers to search for more quantitative relation 
between the fractal dimension of the tumor's contour and 
the efficiency of chemotherapy. Disposing with concrete 
morphometric tumor data, the above inverse geometric 
approach could be used to build the corresponding CA 
model to make predictions of the resistance to chemother- 
apy or provide rough estimates of the therapeutic effi- 
ciency for the respective in vitro grown tumor cell lines. 



Scaling behavior 

Inspired by the fractal and scaling concepts, Bru and 
coworkers [lq . Il9j scrutinized the results of morphome- 
tric experiments on a few tens different in vitro grow- 
ing tumors and cell colonies by the fractal and scaling 
analysis, and suggested that the tumor reveals the basic 
features of Molecular Beam Epitaxy (MBE) universal- 
ity class [l(| characterized by i) a linear growth rate of 
the radius, ii) the proliferating activity at the outer bor- 
der, and iii) diffusion at the tumor surface. Moreover, 
to characterize the tumor growth dynamics, they deter- 
mined the critical exponents |10l. nevertheless the issue 
is still under discussion [2^, |24{. 

Here we outline the typical scaling behavior of the con- 
tours of the growing CA clusters applying the CA rules 
obtained by the above optimization approach for the re- 
spective desired growth rate and fractal dimension. We 
note that no scaling properties have been explicitly in- 
cluded into the optimization procedure. 

A rough interface is usually characterized by the fluctu- 
ations of the contour points rk, n (t), where k refers to the 
A:th point of the nth segment of the arc length I in time t. 

By defining the average r n (t) = (1/N n (t)) Y,k=i r fe,«( i )' 
the local interface width is given by 



w(l,t) 



\ \N n (t) 



N n (t) 

E 

k=i 



\rk, n (t) - r n (t) 



(11) 



where Ni(t), A^t), ... , N n (t), ... are the numbers of 
contour points along the nth angular segment of the arc 
length I. The brackets () mean the average over the dif- 
ferent realizations of the arc length I. 

In general, the growth processes are expected to follow 
the Family- Vicsek ansatz [251 ] . 

u,(i,t) =*"*(*/£(*)), (12) 
with the asymptotic behavior of the scaling function 



u a if u <C 1 
const if ti>l, 



(13) 



where a is the roughness exponent and characterizes the 
stationary regime in which the height-height correlation 
length £0) - t x l z (z is so called dynamic exponent) has 
reached a value larger than The ratio /3 — a/z is 



called the growth exponent and characterizes the short- 
time behavior of the interface. 

Figs. [H] to [TO] demonstrate the typical scaling behavior 
of the growing CA cluster contours. By fitting the slope 
of the w(l,t) in the respective regions (fT3)) of the log- log 
plot we have determined the values of a and (3 exponents, 
giving the value of the exponent z = 0.79. 



0.75 




0.25 



log© 

FIG. 8: Log-log plot of w(l,i). The region with power law 
behavior corresponding to the exponent a shows the slope 
0.5. 



0.75 - 




0.25 



log(t) 

FIG. 9: Log-log plot of w(l,t). The region with power law 
behavior corresponding to the exponent /3 shows the slope 
0.63. 

Performing the scaling analysis of 10 CA clusters grow- 
ing accordingly to the respective CA rules optimized to 
fulfill different combinations of the desired rate constants, 
fractal dimensions and weights w\ we have observed weak 
dependence of the exponent z on the growth rate (Tab.[T|. 
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FIG. 10: Scaling of w(l, t) using the a and z = 0.79 exponents 
as provided by the linear fits (Figs. [8] and [9]). Collapse of the 
curves for large t is obvious. 





A = 0.2 


A = 0.3 


A = 0.4 


a 



z 


0.5 - 0.6 
0.63 - 0.76 
0.79 


0.41 - 0.56 
0.6 - 0.75 
0.63 - 0.74 


0.38 - 0.53 
0.68 - 0.89 
0.55 - 0.6 



TABLE I: The dependence of the scaling exponents values 
on the growth rate of the CA cluster growth. 

The values of the scaling exponents differ from that 
presented in [l9| . Nevertheless, there performed scaling 
analysis of the circular contours, as well as the scaling 
analysis of growin g ci rcular contours in general, is under 
discussion fiH l23l. |24| . 

Conclusion 

During the recent years cancer research has developed 
into a very active scientific field with many concepts 
from other scientific areas, such as statistical physics, 
app lied mathematics, nonlinear science and modeling 
EE M, El, [U, H M, S [H, 11 penetrating into it. 
The ultimate aim of the modeling is, in general, to pro- 
vide correct predictions for the real systems. These tra- 
ditionally come from the deep knowledge of the problem 
(microscopic models, ah initio approaches) or from phe- 
nomenology (e. g. represented by the growth laws writ- 
ten in the form of the first order differential equations). 
The microscopic models of tumor growth are typically 
based on the conceptions such as cell migration, prolifer- 
ation, chemotaxis, apoptosis, etc. However, involving the 
plethora of parameters, these models often face to con- 
ventional problems of complex systems - subjective choice 
of the model parameters and their questionable measure- 
ment, both leading to the loss of predictive power. On 
the other hand, macroscopic models often miss adequacy 



[35l ]. The way out of the obstacle could be based on 
the phenomcnological universalities involving phenomena 
like scale laws, complexity, and nonlinearity |35|. Here 
presented approach can be seen as the instance of that 
way of thinking, benefiting from the demonstrated possi- 
bility to simulate the MTS growth by the respective GA 
optimized CA model. Implementing some of the features 
of the MTS growth into the CA model one can eventually 
observe the feature(s) which has not originally been in- 
cluded into the model and, consequently, try to identify 
the feature (s) in the MTS growth itself. 

To illustrate the future perspective of the above CA 
approach in cancer research we propose the hypothesis 
of the existence of the saturation time dividing the short 
time behavior of the local interface width (fTTj) from its 
universal, time- independent behavior, which is obvious 
in the simulated CA clusters growth (Figs. [8] and [9]). In 
general, during the roughening of a contour the lateral 
correlations develop. At the beginning, different sites on 
the contour are correlated only to the characteristic dis- 
tance, the correlation length (|12p . The local interface 
width increases as a power of time until the correlation 
length reaches the size of the contour, the time is called 
saturation time. After that all the sites become corre- 
lated and the local interface width reaches a (time in- 
dependent) saturation value, which increases as a power 
of the system size. However, if the system, in addition 
to roughening of the contour, increases in size, the ques- 
tion whether and when the saturation is reached (or lost, 
if started as correlated) arises. Consequently, the follow- 
ing questions emerge: Show the different tumors different 
saturation times (and values) of the local interface width? 
Does the behavior of the local interface width relates to 
malignity of tumors? Can it be purposefully affected by 
the chemotherapy? The finding of eventual relation be- 
tween the malignity (or proliferation activity) and the 
saturation time supported by the experimental evidence 
could lead to the development of the new therapeutic ap- 
proach based on the (re)establishing of the correlations. 

The most salient feature of the above approach is that 
neither microscopically based parameters in the tradi- 
tional sense, such as the cell migration, proliferation, 
chemotaxis, apoptosis, etc, are required nor the relations 
between them are a priori built. Both stay hidden in 
the transition rules, optimized to fulfill chosen objectives, 
and acquire statistical (or even symbolic) meaning. For 
example, the proliferation activity can be deduced only 
from the increase (or decrease) of the CA cluster size, 
not from the local changes of the states of individual CA 
sites, as on the here adopted scale it is impossible to dis- 
tinguish between the contributions of cell proliferation, 
migration, apoptosis, etc. On the other hand, some sta- 
tistically modeled mechanisms (such as the response to 
drugs) can be included in a quite straightforward man- 
ner. 
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